Spectral dimension as a probe of the ultraviolet continuum regime of causal 

dynamical triangulations 



Thomas P. Sotiriou, 1,2 Matt Visser, 3 and Silke Weinfurtner, 1 
1 SISSA - ISAS, Via Bonomea 265, 34136, Trieste, Italy and INFN, Sezione di Trieste 

2 DAMTP, CMS, University of Cambridge, Wilberforce Road, Cambridge, CB3 OWA, UK 

3 MSOR, Victoria University of Wellington, PO Box 600, Wellington 6140, New Zealand 

We explore the ultraviolet continuum regime of causal dynamical triangulations, as probed by the 
flow of the spectral dimension. We set up a framework in which one can find continuum theories 
that can in principle fully reproduce the behaviour of the latter in this regime. In particular, we 
show that, in 2 + 1 dimensions, Hofava-Lifshitz gravity can mimic the flow of the spectral dimension 
in causal dynamical triangulations to high accuracy and over a wide range of scales. This seems to 
provide evidence for an important connection between the two theories. 
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The first serious effort to implement discretized ge- 
ometries into the framework of general relativity dates 
back to 1961 and the development of Regge calculus pQ. 
Since then there has been persistent interest in a vari- 
ety of discrete quantum gravity models [3J. One partic- 
ularly interesting variant is causal dynamical triangula- 
tions (CDT) [3] , in which a geometry emerges as the sum 
of all possible triangulations (modulo diffeomorphisms) 
obeying a global time foliation. The sum is evaluated us- 
ing the path-integral formalism, such that every history 
is weighted using a variant of Regge calculus, where the 
edge lengths of the fundamental building blocks, the d+ 1 
dimensional simplices, is kept fixed. The various differ- 
ent histories can be viewed as all possible fluctuations in 
the geometry, and they differ in the number of simplices 
and in the manner the latter are glued together. 

In a discrete model, such as CDT, part of the chal- 
lenge is to find suitable probes for the continuum limit. 
The latter has been investigated numerically, working at 
a fixed volume (fixed number of simplices) and using a 
Monte Carlo algorithm [H [5]. The spectral dimension 
has been proposed as a probe in CDT and has a longer 
history in lattice quantum gravity, see Refs. [El [7] and ref- 
erences therein. This can be thought of as the effective 
dimension as probed by an appropriately defined (ficti- 
tious) diffusion process (random walk) . More concretely, 
a diffusion process from point x to point x' in (fictitious) 
diffusion time s is characterized by the probability den- 
sity p(x, x',s). The average return probability P(s) is 
defined as /?(x, x, s) averaged over all points in space. 
The spectral dimension is then defined as 



d, = -2 



d In P(s) 
dln(s) 



(1) 



In CDT the diffusion is a discretized, stochastic pro- 
cess. The spectral dimension is not constant, but instead 
changes with s, and consequently with the length scale. 

Clearly, the concept of the spectral dimension is not 
limited to CDT or discrete theories. One could consider 
a diffusion process on a manifold, as was proposed in 



Ref. [7]. Then the spectral dimension is defined same 
as above and p(x, x',s) is determined by the diffusion 
equation 



dp(x,x',s) 
ds 



+ Dp(x,x',s) = 0. 



(2) 



The choice of the differential operator D corresponds to 
the "type" of diffusion process being considered. 

For example, when D is the 3-dimensional Laplacian 
and s is identified with real time Eq. ^ becomes the heat 
equation. Instead, when a diffusion process is to be used 
as a probe of geometry (or kinematics), s becomes a ficti- 
tious diffusion time and x represents a point in spacetime 
(in analogy with what was mentioned above for CDT). 
A natural choice for D then seems to be the propagator 
of perturbations in spacetime. For instance, one could 
have D — g^VaV^, where g a p is the (Lorentzian) met- 
ric on a manifold (Greek indices run from to d, which 
denotes the number of spatial dimensions) and we have 
performed a Wick rotation (t —> —it). In this case, for 
large s the corresponding spectral dimension will probe 
the geometry associated with g a fj. For small s, given that 
spacetime is flat in a sufficiently small neighborhood of 
each point, the spectral dimension will actually probe the 
kinematics associated with D [8] . So, when D reduces to 
the flat d'Alcmbertian at small scales, then d 3 — > d + 1. 
Note that if D is a more complicated differential oper- 
ator at small scales {e.g. due to ultraviolet corrections) 
this will be encoded in the small s behaviour of d s . 

Given the above, it is very tempting to use the spectral 
dimension as a probe of the continuum limit of discrete 
theories, and CDT in particular, or as a potential link 
to continuum theories. In fact, in Ref. [5] the large s 
behaviour of d s was matched to the outcome of a dif- 
fusion process in a curved manifold with the geometry 
of a stretched sphere. However, the small s (ultravio- 
let) behaviour of d s in CDT is, to date, not well un- 
derstood. This regime is believed to encode the influ- 
ence of quantum corrections and is expected to provide 
the most prominent hint towards the effective field the- 
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ory arising from CDT. Our goal is to demonstrate that 
one can indeed extract considerable information regard- 
ing the ultraviolet continuum regime of CDT from the 
small s behaviour of d s , and make a first crucial step 
towards identifying the characteristics of an (effective) 
continuum theory that could reproduce this behaviour. 

The key idea is that d s (s) can be used to determine D, 
which in turn characterizes (to some extent) a continuum 
theory. In Ref. jB] we argue that in principle one can 
indeed determine the dispersion relation associated with 
D when d s (s) is known in closed form. However, when 
the latter is known only in tabulated form, as is the case 
for CDT, one should instead rely on the techique of non- 
linear regression to reconstruct the dispersion relation 
from d s (s). This requires some educated guess for the 
general form of the dispersion relation. 

Consequently, what one can do in practice is to choose 
a continuum theory of gravity that shares some funda- 
mental characteristic (s) with CDT and check if it can re- 
produce the behaviour of d s (s). Clearly, this would not 
prove that the theory in question is the continuum limit 
of CDT: the spectral dimension does not carry all the in- 
formation about the theory (as we will see in more detail 
below). Additionally, due to the differences in the defini- 
tions of the spectral dimension in the discretium and in 
the continuum there is still some ambiguity on the nature 
of the linkage it provides between discrete and continuum 
theories. Nonetheless, even with these caveats in mind, 
identifying a continuum theory that could fully repro- 
duce the behaviour of the spectral dimension for small 
s in CDT certainly seems to be a crucial step in under- 
standing its ultraviolet continuum regime. 

A suggestion for the candidate continuum theory has 
been made in Ref. [7 . It was noticed there that the 
preferred causal structure of CDT, imposed by a pre- 
ferred foliation by slices of constant time, is reminiscent 
of the symmetries of Hofava-Lifshitz (HL) gravity [S]. 
(See Ref. [TU] for a brief review including viability con- 
straints.) Indeed, the latter is a theory with a preferred 
spacelike foliation, described by a scalar field. The ex- 
istence of the preferred foliation allows for higher order 
spatial derivatives in the theory without having higher or- 
der time derivatives. This leads to significantly improved 
ultraviolet behaviour which actually renders the theory 
power-counting renormalizable [9l lllj , at the expense of 
giving up Lorentz invariance. 

The action of HL gravity is, in the preferred foliation, 



S = 



M 2 pl 



J d d xdtN^/g{K ij K i:j - XK 2 



(3) 



where K io = - ViNj - VjNi) /(2N) is the extrin- 
sic curvature, N is the lapse, Nj the shift and gij is the 
induced metric on the spacelike hypersurfaces (Latin in- 
dices run from 1 to d). M p \ is the Planck mass and A is a 
dimensionless running coupling. V is the part of the La- 
grangian which contains only spatial derivatives. Power- 



counting renormalizability requires that it includes terms 
with at least 2d derivatives. Generally, V should include 
all terms compatible with the symmetries of the theory 
[H [T3]. (For A = 1, V = R, where R is the Ricci scalar 
of g^, the theory reduces to general relativity.) 

We will focus here on 2 + 1 dimensions (d = 2) mainly 
because, due to computational limitations, lower dimen- 
sional simulations are expected to produce more accurate 
data sets for our purposes. For d = 2 we have [T3] 

V = £R + T]a 2 + 9l R 2 +g 2 V 2 R + g 3 a 4 +g 4 Ra 2 

+g b a 2 {Wa 3 ) + g 6 (V l a t ) 2 + ff7 (V ^(VV), (4) 

where = dAnN, a 2 = a^a* and r\ and £ are dimension- 
less coupling, while the gi couplings have dimensions of 
an inverse mass squared. We have neglected the cosmo- 
logical constant as it is irrelevant for our purposes. 

We are interested in small length scales where curva- 
ture effects are negligible and quantum effects are ex- 
pected to be important. So, for the diffusion process 
through which we define d s for HL gravity it is sufficient 
to use linearized propagators around flat space for D. As 
we have shown in Ref. |14j , in 2 + 1 dimensional HL grav- 
ity, the foliation-defining scalar with dispersion relation 



Pi(k 2 
P 2 (k 2 
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1+Dk 2 



is the only propagating mode. Here 
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(The dispersion relation for the scalar mode is 3 + I di- 
mensional HL gravity is qualitatively the same |13j.) 

In Ref. [5] we discuss in detail how one can define the 
spectral dimension associated with a general dispersion 



relation of the form 



f(k ). In this case the formal 



solution of Eq. ^ after Wick rotation is 



p(x,x', 



d 2 kdu 
(2tt) 3 ' 



i(k-(x-x')) „-s(w 2 +/(fc 2 )) 



A straightforward calculation yields 



4 = 1 + 2s 



/ /(fc 2 ^-^- 8 ^ 2 ) dfc 

j k d-l e -3 f(V) dk 



(7) 



(8) 



This formula is directly applicable to the dispersion re- 
lation in Eq. (pb. 

We have now laid out all of the technical tools that 
allow us to fit the spectral dimension for CDT in the 
ultraviolet continuum regime in 2 + I dimensions using 
the simulations of Ref. [5]. These simulations were re- 
peated for 6 different values of the number of simplices 
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N s = [40,50,70,100,140,200] x 10 3 . For each indepen- 
dent simulation N s was held fixed, and one sums over 
1000 different histories (geometries). On each history the 
(fictitious) diffusion process was implemented and used 
to calculate the corresponding spectral dimension d s . 

In Ref. |7j the spectral dimensions of CDT and HL 
gravity were compared at the two limits of the ultravio- 
let continuum regime and they were found to be in good 
agreement. Both theories seem to predict that at small 
scales d s flows to 2, whereas at length scales large enough 
for the quantum correction to become largely subdomi- 
nant, but small enough for curvature correction to be 
unimportant, d s flows to 3, the number of topological di- 
mensions. Encouraging as it may be, this result is far 
from being conclusive, given that there is an infinity of 
curves joining two points. 

Our goal is far more ambitious: we numerically eval- 
uate the spectral dimensions in Eq. {[sj) for a dispersion 
relation of the kind ([5]), and fit the CDT data using the 
method of non-linear regression. This will allow us to 
obtain the parameters of the dispersion relation that pro- 
vide the best fit for the whole ultraviolet continuum re- 
gion, and access how good such a fit can be. This will 
provide a conclusive answer to whether HL gravity can 
reproduce the behaviour of d s for CDT in this regime. 

This procedure will not uniquely determine all of the 
parameters of HL gravity. First of all, it is already clear 
from Eqs. ([6| that g$, g\ and g§ do not enter into the 
linearized propagator and, therefore, they cannot be de- 
termined without taking into account curvature correc- 
tions. Secondly, there are 7 couplings that do enter in 
the propagator, but only in 4 combinations. This implies 
that any observable that probes the dispersion relation 
will give us some but not full insight into the fundamental 
theory. Finally, without loss of generality, we can choose 
to work in units where the infrared light speed is one, 
i.e. A = 1. This amounts to the rescaling u —> u)j\fk. 

In Fig. 1(a) we plot the spectral dimension d s as re- 



sulting from the CDT simulation with the largest number 
of simplices. As illustrated with the differently shaded 
areas, the diffusion process can be divided into three 
physically different regions. For very small s £ [1 , 15] 
the spectral dimension shows discrete behaviour as ex- 
pected. For very large s £ [614, 29999], the behaviour 
of the spectral dimension is expected to be dominated 
by curvature effects. This is the region were the flow of 
d s is compatible with a diffusion process on a continu- 
ous manifold with the geometry of a stretched sphere [5] . 
The intermediate region, where a continuous behaviour 
emerges but quantum corrections are pre-dominant, is 
the one we are interested in. Our model-fit for this re- 
gion (s £ [15, 614]) is represented by a red solid line. 
Fig. 1(b) zooms into this region. 

We have chosen to present a fit for the simulation with 
the larger number of simplices, because the higher the 
number of simplices, the larger the volume of the uni- 



verse will be in the simulations. If the volume is not large 
enough, curvature effects can kick in before quantum 
effects become subdominant. There might then be no 
regime where both effects can be neglected, so no regime 
where the diffusion effectively takes place in flat space 
and the spectral dimension flows to 3. This undesirable 
behaviour is indeed observed in some of the simulations 
with fewer simplices. 



In Fig. 1(a) we see that for the 200, 000 simplices simu- 
lation d s starts at around 2 in the ultraviolet continuum 
regime, it reaches 3 in a regime where quantum effects 
are apparently still non-negligible and push d s slightly 
over 3. Eventually, d s reaches 3 for a second time and 
then drops to lower values under the influence of cur- 
vature corrections. Remarkably, our model can repro- 
duce this behaviour very well in all the region where cur- 
vature corrections are negligible, including the part of 
the curve that exceeds 3. The corresponding values for 
the parameters of the model, in units where A = 1, are 
B = -1.18 ± 0.22, C = 344.47 ± 3.11, D = 10.08 ± 0.30, 
which are compatible with technically natural values for 
the various couplings of HL gravity. 

As an indication of how good our fit is, we present 



in Fig. 1(c) the individual residuals (difference between 



data and best fit) for 3 simulations involving different 
number of simplices: N s = [50, 100, 200] x 10 3 . For all 3 
cases the residuals grow for very small values of s. We 
attribute this to the fact that at some sufficiently small 
s the discrete nature of CDT becomes important. Note 
that, even though we do not attempt to fit the lightly- 



shaded area of Fig. 1(a) where the data points violently 
oscillate, the absence of obvious oscillations in d s does 
not imply the complete absence of discreteness effects. 

Moving to intermediate values of s, a somewhat wor- 
risome feature of the residuals is that they are not ran- 
domly scattered, but instead they exhibit some oscilla- 
tory pattern. This indicates that the fit might be missing 
a systematic effect. Nevertheless, the fact that the ampli- 
tude of the oscillations clearly reduces significantly as one 
moves to larger number of simplices is very encouraging, 
as it indicates that this effect is likely to become negligi- 
ble as the size of the simulations increases. Additionally, 
this unknown effect might well be the subdominant con- 
tribution of the aforementioned discreteness effects. 

The absolute magnitude of the residuals for all sim- 
ulations becomes large again for large enough values of 
s. This is expected and signals the scale at which cur- 
vature corrections, which we neglect, become important. 
Consequently, the fact that the residuals start growing 
at a value of s that is lower that the value for which 
d s = 3 for the second time (or even for the first time 
for lower numbers of simplices) seems to indicate that, 
even for the largest number of simplices, the universe is 
not large enough to allow one to comfortably define a 
patch of spacetime that is large enough for quantum ef- 
fects to be negligible, and at the same time small enough 
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FIG. 1: (a) The spectral dimension as a function of (fictitious) diffusion time s (semilog scale). Crosses represent data points 
for CDT for simulation with 200, 000 simplices, whereas the red solid line represent the fit based of HL gravity. The region 
where discretium effects are important is lightly shaded. The region where curvature effects dominate is darkly shaded. The 
intermediate is the ultraviolet continuum region, (b) Zoom into the latter region, (c) Individual residuals for the model fits 
for 50, 000, 100, 000 and 200, 000 simplices. The residuals grow as expected at very small and very large scales under the 
influence of discretium and curvature effects respectively. The oscillatory pattern at intermediate scales seems to indicate some 
systematic deviation from the model, which however, tends to disappear as the number of simplices is increased. For 200, 000 
simplices the fractional residual remains below 0.5% almost throughout the entire ultraviolet continuum region. 



for curvature corrections to be largely subdominant. 

It should be pointed out that we have neglected the 
running of couplings, even though d s is not necessarily 
a universal observable throughout the region considered. 
We are not fitting the deep UV, and so it is not unreason- 
able to assume that the running is not very significant. 
Our results a posteriori justify this expectation. 

We have attempted to fit the CDT data with other 
types of dispersion relations, either inspired by lattice 
field theory, or just simple polynomials, which typically 
arise in projectable HL gravity where the lapse function 
is taken to be space independent [5J[T2]. However, the 
independent residual test largely favours rational disper- 
sion relations, as in Eq. ^ [IS]. The flow of d s as defined 
in the discretium in CDT appears to match that of the 
foliation-defining scalar in HL gravity in the ultraviolet 
continuum regime, and not that of a test scalar field. 

To summarize, we have shown that in 2 + 1 dimen- 
sions HL gravity provides a very good fit (which is ex- 
pected to improve as simulations include a larger num- 
ber of simplices) to the flow of the spectral dimension in 
CDT throughout the ultraviolet continuum regime. Even 
though this is by no means enough to argue that the for- 
mer is the continuum limit of the latter, it does provide 
evidence for an important connection between the two 
theories. It also demonstrates beyond doubt, that the 
spectral dimension is a powerful tool for relating discrete 
and continuum theories and for gaining insight into their 
behaviour. 
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